Mass dependence of light nucleus production in ultrarelativistic heavy ion collisions 
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Light nuclei can be produced in the central reaction zone 
via coalescence in relativistic heavy ion collisions. E864 at 
BNL has measured the production of ten light nuclei with 
nuclear number of A = 1 to A = 7 at rapidity y ~ 1.9 and 
Pt I A < 300MeV/c. Data were taken with a Au beam of mo- 
mentum of 11.5 A GeV/c on a Pb or Pt target with different 
experimental settings. The invariant yields show a striking 
exponential dependence on nuclear number with a penalty 
factor of about 50 per additional nucleon. Detailed analysis 
reveals that the production may depend on the spin factor of 
the nucleus and the nuclear binding energy as well. 

Relativistic heavy ion collisions may create high energy 
density and high baryon density in the reaction zone. 
Therefore, they are considered to be the best laboratory 
environment to create novel objects or states such as the 
Quark-Gluon Plasma @. At the AGS energies (y/s ~ 4.8 
AGeV), the baryon stopping is high, which means that 
the baryon density near center-of-mass rapidity is quite 
large compared to the much higher beam energy at the 
SPS (Vs ~ 17 AGeV) or at the RHIC {yfs ~ 200AGeV") 
Light nuclei can be produced by the coalescence of 
created or stopped nucleons j|. Since the probability of 
coalescence of a particular nuclear system (d, 3 He, etc.) 
depends on the properties of the hadronic system formed 
as a result of the collision, the study of the coalescence 
process is useful in elucidating those properties. For ex- 
ample, in a coalescence model, the coalescence probabil- 
ity depends on the temperature, baryon chemical poten- 
tial (essentially the baryon density), and the size of the 
system, as well as the statistical weight of the coalesced 
nucleus fl6[] . A thermal model |l5) does not use details 



of cluster creation, however, the dependence of the pro- 
duction on the system configuration is very similar to 
the coalescence model. The data presented in this paper 
shows evidence that the probability may also depend on 
the binding energy of the coalesced nucleus. Systematic 
study of the production of light nuclei is limited by their 
low production rates in relativistic heavy ion collisions. 
E864 is the only experiment which is able to measure the 
production of charged nuclei with A > 4 produced in the 
central reaction zone. However, the nature of the coales- 
cence process tells us that the higher the baryon number, 
the more sensitive the production rate is to the system's 
configuration. 

In this letter, the production of protons, neutrons, 
deuterons, 3 He, tritons, 4 He, 6 He, 6 Li, 7 Li and 7 Be 
around rapidity y ~ 1.9 and transverse momentum of 
p T /A < 300MeV7c in 10% most central Au+Pt(Pb) col- 
lisions measured by E864 Q at BNL is presented. These 
measurements have significant impact on the strange 
quark matter Q search by several experiments in rela- 
tivistic heavy ion collisions. They also provide informa- 
tion about the thermal equilibrium of the system and the 
detailed process of coalescence. 

E864 H at Brookhaven National Laboratory is an open 
geometry, high data rate apparatus designed to search 
for novel objects such as strange quark matter. It is 
a fixed-target experiment at the AGS with incident Au 
beam momentum of W.bGeV/ c per nucleon. We can have 
two independent mass measurements with the E864 ap- 
paratus :the tracking system and the hadronic calorime- 
ter. The requirement that the two measurements agree 
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provides excellent particle identification and background 
rejection. The tracking system has two dipole analyz- 
ing magnets (Ml and M2) followed by three hodoscope 
planes (HI, H2 and H3). There are two straw stations 
(S2 and S3), each with three close-packed double planes. 
Each scintillating hodoscope plane has 206 vertical slats, 
and there are 960 and 1920 straw tubes in each S2 and S3 
straw plane respectively. The tracking system measures 
momentum, charge and velocity (/3) of charged particles 
with a mass resolution of few percents in the region of in- 
terest. Charge misidentification is less than 1 in 10 8 due 
to three redundant measurements of energy loss in the 
hodoscopes. There is an additional mass measurement 
from the hadronic calorimeter || which has good energy 
(AE/y/E = 0.344/v 7 ^ + 0.035) and time (a t ^ 400ps) 
resolutions. It is made of 754 towers of scintillating fiber- 
embedded lead. The containment of a hadronic shower is 
about 50% in the peak tower and about 90% in the 3x3 
array around the peak tower. This allows sufficient qual- 
ity requirements of the shower cluster and provides for 
good neutron measurements. A vacuum tank along the 
beam line reduces the background from beam particles 
interacting with air. The total length of the apparatus is 
about 28 meters. 

The trigger consists of good beam definition, a multi- 
plicity requirement [JtJ and an optional level II high mass 
trigger - the Late-Energy Trigger (LET) f§. The LET 
utilizes the energy and time-of-flight measurements from 
the calorimeter. Only the 10% most central collisions 
were collected for the data presented here. The data 
are from four different experimental settings optimized 
for the different physics topics. Proton, deuteron, 3 He 
and triton measurements are from the 1995 run with the 
"+.45T" field setting on Ml and M2 magnets @ and 
with a Pb target of 5% interaction length for Au nuclei. 
Neutron data are from the 1995 run with the highest field 
setting of " + 1.5T" (5% Pb target) (l3| while the rest of 
the data are from the 1996 run at the " + 1.5T" field set- 
ting with the high mass level II trigger and Pt target 
of 60% interaction length (physical thickness of 1.5cm) 
mainly for the strange quark matter search with the 
exception that the Be and Li measurements were from 
the combination of "-.75T" and "+1.5T" field settings 
taken in the 1996 run. Due to the open geometry, there 
is sufficient overlap from different settings for consistency 
checks @-[l3). 

Because of E864's large acceptance, we are able to 
divide each measurement into several rapidity (y) and 
transverse momentum (pt) bins. At y ~ 1.9 andpy/A < 
300MeV/c, the acceptance is about 20% for most of the 
light nuclei. It is in this momentum range where we can 
compare yields of the different species. In other bins, we 
do not have acceptance for all the light nuclei detected. 
Interesting and more detailed aspects of the rapidity and 
transverse momentum distribution of the produced light 
nuclei will be presented in other papers j|. 

For charged nuclei with A < 3, the production rates 



are relatively high so that the LET trigger is not required 
and the analysis does not use the hadronic calorimeter. 
The calorimeter is the only significant detector for the 
neutron measurement where the tracking detectors serve 
as tools to calibrate the system and to reject the charged 
particles which deposit energy in the calorimeter. For 
high mass nuclei (A > 4) with low production rate, the 
LET was required to select those events with high mass 
candidates. The LET rejects those events without any 
high mass candidates and achieves a rejection factor of 
about 70 while maintaining good efficiency for high mass 
states (~ 50% for 4 He and >85% for A > 5). 

Invariant yields are calculated and presented in terms 
of d 2 NJ (2irpTdpT<iy) in units of GeV~ 2 c 2 as shown in 
Table §. The number of particles (N) are taken from the 
mass distributions in each rapidity and px bin with back- 
ground subtracted. The background level is about 10% 
or less in most of the bins. Acceptance and efficiency 
corrections using either simulation or data analysis are 
applied with the associated systematic and statistical er- 
rors. The errors of the measurements of proton, neutron, 
deuteron and 3 He are on the order of 10%. The rest of 
the measurements have larger errors varying from 20% 
(a, e He) to 45% ( 6 Li, 7 Li, 7 Be) due to uncertainty in 
the LET trigger efficiency (~ 10%), low statistics (up to 
40%) and high background level from lower mass nuclei 
(up to 20%). Target absorption of the produced parti- 
cles is negligible for the 5% Pb target and about 8% for 
the 60% Pt target. The total systematic error for these 
measurements varies from a few percent to about 25%. 

Figure [l] shows the invariant yield as a function of nu- 
clear number A for stable or metastable particle with 
A = 1 to A = 7. The rapidity binning is Ay = 0.2 
with the exception of Ay = 0.6 for 1 Li. px bin- 
ning is Ap T = lOOMeV (100MeV-200MeV) for A=l 
to 3, 250Me^ (500MeV-750MeV) for 4 He, 500Me^ 
(500MeV-lGeV) for 6 He and 6 Li, and 2GeV (0-2GeV, 
acceptance peaks at 750MeV) for 7 Li and 7 Be. These 
bin sizes keep pr/A < 300AfeV. The invariant yields 
span almost ten orders of magnitudes with striking ex- 
ponential behavior. 

A fit to the A dependence of the invariant yields in 
this rapidity bin results in a penalty factor of about 50 
(48 ± 3, x 2 /8 = 4.9) for each additional nucleon to the 
nuclear cluster. This penalty factor is much higher than 
the penalty factor in the system with lower beam energy 
at the BEVALAC The consequence is that it is 

much harder to form high mass objects by coalescence, 
such as strange quark matter. 

In a statistical approach to the formation of light nu- 
clei, the yield is proportional to the spin factor (2J+1) 
SO "III- One can go one step further to analyze the 
deviations of the invariant yields from the exponential 
behavior (figure §.A). The measured ratios of proton to 
neutron, 3 He to triton, 6 Li to 6 He and their correspond- 
ing spin factors strongly indicate that the production rate 
is proportional to the spin factor (2J+1) of the produced 
particle ||,|ll| |l3|. Therefore, it is probably reasonable to 
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include this spin factor in the production rate, which is 
consistent with most of the models [^|Jl^-|l^1. Figure [|.B 
shows the spin corrected deviations (from exponential be- 
havior), which still have significant deviations. We note, 
however, that the ratios for A — 6 states with spin fac- 
tors differing by a factor of 3 are brought into agreement 
with each other by the spin factor correction ]g|Jl^|. The 
spin correction factor is taken as (2 J + l)/(2 x| + 1) 
with a normalization to the nucleon spin factor of 2. 

When the deviation after spin factor correction is plot- 
ted as a function of binding energy per baryon as shown 
in Figure ||, we can fit the dependence with an exponen- 
tial function of an inverse slope of T s = 5.9 ± l.lMeV 
when the A dependence of (2 J + l)/2 x 26/48 A ~ 1 is ap- 
plied. The small difference (n to p ratio of about 1.2±0.1) 
between the abundances of neutron and proton at freeze- 
out J|,|ll],[l3| is corrected for in the analysis shown in 
Figure If the total binding energy instead of the 
binding energy per baryon is used, the inverse slope is 
about T s ~ 36MeU. 

From the measurements of light nuclei production near 
midrapidity (at y ~ 1.9, where ycM — 1-6) with low 
transverse momentum, a penalty factor of about 50 for 
each additional nucleon is found in the invariant yield. 
Although the total production rate comes from the in- 
tegration of the whole phase space which might differ 
from the measured invariant yield (due to flow, etc.), we 
do not expect the ratios of different particle species be- 
tween total production and the rapidity and pt range we 
cover to be different by orders of magnitude. In fact, if 
we use the parametrization of the correction factor from 
iPI , the penalty factor is estimated to be between 39 
and 72. Therefore, the penalty factor can be applied to 
estimate the production rates of heavy clusters. Due to 
the small elastic structure function of light nuclei at large 
momentum transfer |19|, the possibility of light nuclei at 
1. < y < 2.2 coming from projectile or target is small 
compared to the observed production rates. 

Naive comparison with theoretical estimates of chem- 
ical potential and temperature at freeze-out from low 
mass hadronic spectra [jl5| will show us whether the 
light nuclei might have different chemical potential and 
temperature at hadronic freeze-out. If we use the 
chemical potential of /ijy = 540MeV and tempera- 
ture of T = 120MeV [jl5|], we get a penalty factor of 
exp [(mjv — (J-n)/T] ~ 28. A kinetic freeze-out and ra- 
dial flow analysis shows a higher penalty factor of about 
75 with n N = 536AfeU and T ~ 93MeV @. 

The dependence of production rate on binding energy 
per baryon shows the sensitivity of our data and can not 
be explained by the coalescence model or the thermal 
model with the simple exp [-B/T] where T ~ lOOMeV 
and B is the total binding energy [jl7|. One possible 
explanation is that the production rate might depend on 
the size of the produced object, and therefore depend on 
its binding energy. It is also possible that there exists 
some subtle final state interaction which depends on the 
size or the binding energy. For example, collisions with 



sufficient energy to break up the nuclei can occur down 
to surviving temperatures which are comparable with the 
binding energy [ p0[ . 

In summary, E864 measures the invariant yields of 
light nuclei production of protons, neutrons, deuterons, 
3 He, triton, 4 He, 6 Li, 6 He, 7 Li and 7 Be near ycM and 
Pt — 0. A striking exponential behavior of light nucleus 
production as a function of nuclear number is found. The 
penalty factor, extracted from the light nuclei production 
rates, for an additional nucleon in the nuclear cluster is 
about 50. Detailed analysis reveals that the invariant 
yield may also depend on the spin factor and the binding 
energy of the produced nucleus. 

We gratefully acknowledge the excellent support of the 
AGS staff. This work was supported in part by grants 
from the U.S. Department of Energy's High Energy and 
Nuclear Physics Divisions, and the U.S. National Science 
Foundation. 

* Present address: (11) 

f Present Address: Anderson Consulting, Hartford, CT 
X Present address: Univ. of Denver, Denver CO 80208 
§ Deceased. 

|| Present address: McKinsey & Co., New York, NY 10022 

% Present address: Department of Radiation Oncology, Med- 
ical College of Virginia, Richmond VA 23298 

t] Present address: University of Tennessee, Knoxville TN 
37996 

j) Present address: Institut de Physique Nucleaire, 91406 Or- 
say Cedex, France 

t Present Address: Institute for Defense Analysis, Alexandria 
VA 22311 

ft Present Address: MIT Lincoln Laboratory, Lexington MA 
02420-9185 

Present address: (2) 



[1] J. W. Harris and B. Muller, Ann. Rev. Nucl. and Part. 

Sci. 46(1996)71. 
[2] T. Wienold et al., Nucl. Phys. A 610 (1996) 76c. 
[3] J. L. Nagle et al., Phys. Rev. C 53(1996) 367, references 

therein 

[4] R. L. Jaffe, Phys. Rev. Lett. 38(1977)617E. 
[5] T. A. Armstrong et al., Nucl. Inst. Meth. A406(1998)227. 
[6] T. A. Armstrong et al., accepted for publication in Nucl. 
Inst. Meth. A 

[7] P. Haridas et al., Nucl. Inst. Meth. A385(1997) 412. 

[8] J. C. Hill et al., Nucl. Inst. Meth. A421(1999)431. 

[9] Manuscript in preparation, see also llCj — P_3|] . 
[10] J. K. Pope, et al., pll9, Proceedings of Heavy Ion Physics 
at the AGS (HIPAGS'96, edited by C. A. Pruneau, et al, 
Published by Wayne State Univeristy); J. K. Pope, PhD 
Thesis, Yale University 1997. 



3 



[11] 



[12] 



[13] 



[14] 
[15] 



[16] 

[17] 

[18] 
[19] 

[20] 



N. George, et al., Proceedings of 15th Winter Work- 
shop on Nuclear Dynamics, Park City, UT, January 13, 
1999(edt. W. Bauer, et al., Kluwer Academic/plenum 
Publishers); N. George, PhD Thesis, Yale University 
1999. 

Z. Xu, et al., J.Phy s. G25 (1999) 403-410; E. Finch, et 
al., [nucl-ex/9908019[ Z. Xu, PhD Thesis, Yale University 
1999. 

T. A. Armstrong et al., nucl-ex/9909001 , accepted for 
publication in Phys. Rev. C; E. Finch, PhD Thesis, Yale 
University 1999. 

A. Baltz et al., Phys. Lett. B325(1994) 7. 
P. Braun-Munzinger and J. Stachel, J. Phys. G: Nucl. 
Part. Phys. 21 (1995) L17; P. Braun-Munzinger et al., 
Phys. Lett. B344(1995) 43. 

R. Sc heibl et al., Phys. Rev. C 59(1999) 1585; H. Dobler, 
et al. 



nucl-th/990401£ 



For example, S. Das Gupta and A. Z. Mekjian, Phys. 
Rep. 72 (1981) 131. 

S. Nagamiya et al., Phys. Rev. C 24 (1981) 971. 

L. C. Alexa et al., Phys. Rev. Lett. 82(1999) 1374; D. 

Abbott et al., Phys. Rev. Lett. 82(1999) 1379. 

V. Serfling, et al., Phys. Rev. Lett. 80(1998) 3928. 




FIG. 1. Invariant yield as a function of nuclear number A. 
The bin size in rapidity is 0.2 (0.6 for 7 Li). p T /A < 300MeV. 
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FIG. 2. Ratio of the invariant yield to the exponential func- 
tion vs. nuclear number A with and without spin factor. The 
top panel is the ratio and the bottom panel is the ratio with 
spin correction. See text for details. 
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Binding Energy per nucleon (MeV) 

FIG. 3. Ratio of the spin-isospin-factor-corrected invariant 
yield to the exponential function exhibits exponential depen- 
dence on the binding energy per nucleon. The isospin abun- 
dance is (n/p)=1.2. Alphas have the largest binding energy 
per nucleon and therefore have the highest production rate de- 
viation from the global exponential behavior. Neutrons and 
protons are not in this plot. Binding energies are calculated 
according to the mass difference between the nucleus and its 
constituent nuclcons. See text for details. 
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TABLE I. Invariant yields of light nuclei production at y — 1.9 and pr/A < ZQQMeV/c. The bin size in y is Ay — 0.2. 
Invariant yields are calculated in terms of d 2 N/ (2-irpTdpTdy) in units of GeV~ 2 c 2 . See details in text. 



Species 


P 


n 


d 


t 


3 He 


Yield 




29.0 


25.8 


0.567 


1.04 x 10~ 2 


9.65 x 10~ a 


((GeV/c)- 


2 ) 


±3.2 


±1.7 


±0.034 


±6.6 x 10~ 4 


±3.3 x 10 -3 


Species 


4 He 


b He 


«Li 


7 Li 


1 Be 


Yield 




2.55 x 10 -4 


5.2 x 10 -8 


2.1 x 10"' 


0.92 x 10 -8 


1.3 x 10~ 8 


((GeV/c)- 


2 ) 


±6.7 x 10~ 5 


±1.2 x 10~ 8 


±7.9 x 10~ 8 


±4.0 x 10~ 9 


(±3.5-2.5) x 10" 9 
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